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We study strategies for estabUshing long-distance entanglement in quantum networks. Specifically, 
we consider networks consisting of regular lattices of nodes, in which the nearest neighbors share a 
pure, but non-maximally entangled pair of qubits. We look for strategies that use local operations 
and classical communication. We compare the classical entanglement percolation protocol, in which 
every network connection is converted with a certain probability to a singlet, with protocols in 
which classical entanglement percolation is preceded by measurements designed to transform the 
lattice structure in a way that enhances entanglement percolation. We analyze five examples of such 
comparisons between protocols and point out certain rules and regularities in their performance as 
a function of degree of entanglement and choice of operations. 

PACS numbers: 03.67.-a, 03.67.Bg, 64.60.ah 



INTRODUCTION 



Entanglement is the property of states of multipartite 
quantum systems that is the most important resource 
for quantum information processing jlj . For this reason, 
one of the most important tasks of quantum information 
science is to establish entanglement at long distances in 
quantum networks, and to optimize final entanglement 
and probability of success. 

Quantum networks 2, 3] play a key role in quantum 
information processing. Here we limit our attention to 
those networks in which quantum states can be prepared 
initially and shared. That is, entanglement can be gen- 
erated between neighboring or, at least, not-too-widely- 
separated nodes (or stations). There are two instances in 
which the above-mentioned tasks become obviously rele- 
vant. On one hand, one can consider macroscopic quan- 
tum communication networks, such as cryptographic net- 
works, or more generally quantum communication nets 
[J, [a, Q, or distributed quantum computation [3| in- 
volving arbitrary nodes of the network. The second in- 
stance concerns microscopic or mesoscopic networks that 
could constitute architectures of quantum computers (cf. 
Ref.i). 

Despite enormous progress in experimental techniques 
(cf. Ref. la and references therein), it is in principle a 
very hard task to establish entanglement at large dis- 
tances due to decoherence and attenuation effects. Two 
remedies for this problem have been proposed: 



• Quantum repeaters. This concept has been de- 
veloped for ID quantum communication chains 
|9l. KLOl. [ill [l3l . Although the simple entanglement 
swapping [l3| procedure can lead to quantum com- 
munication at large distances (see Fig. [T]), for im- 
perfect resources, the performance of such commu- 
nication chains decays exponentially with the dis- 
tance (i.e. the number of repeaters). However, one 
can use more sophisticated quantum repeater pro- 
tocols, which use purification and swapping meth- 
ods that lead to polynomial decay only. 

• Entanglement percolation. Recently, our collabora- 
tors, together with two of us, proposed using net- 
works in which properties of the connectivity of 
the network enable the establishment of, and de- 
termine the probability of, entanglement on large 
distances. In Refs. [ij and [la we considered in 
particular pure-state networks on regular lattices, 
where the nearest-neighbor (NN) nodes share a 
non-maximally entangled pair of qubits, or more 
generally qudits (an entangled bond) . We searched 
for local operations and classical communication 
(LOCC) protocols that lead to establishment of en- 
tanglement between remote nodes of the network. 

At present, only a handful of studies that explicitly 
apply percolation theory to problems in quantum infor- 
mation have appeared in the literature. These articles, 
comprising Refs. uA and [la as well as applications to clus- 
ter states, are reviewed in Ref. [l^. 

Our results leading to the present work can be sum- 
marized as follows: 
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For ID chains we proved that even optimal LOCC 
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FIG. 1: Entanglement swapping and ID repeater. Circles are 
qubits. Heavy lines represent pure non-maximally-entangled 
states. Loops represent the entanglement swapping measure- 
ment on a pair of qubits. a) Two states a and f3 of the form 
specified in ([T]). After the operation, qubits a and d may be 
in an entangled, mixed state, b) swapping as the first step in 
a repeater, c) a ID chain of repeaters. 



strategies only allow the establishment of entangle- 
ment between distant nodes that decays exponen- 
tially with both the distance, as well as the quality 
of entanglement of NN bonds. 

• In 2D and higher dimensions, the possibilities for 
protocols are greatly expanded. The most straight- 
forward, naive protocol — the one we use as a base- 
line to evaluate other protocols — is the one we 
term classical entanglement percolation (CEP), al- 
though it does involve obviously some quantum op- 
erations. This protocol begins with converting (us- 
ing LOCC) each of the entangled bonds into a sin- 
glet (i.e. maximally entangled state) with proba- 
bility p (the so-called singlet conversion probabil- 
ity (SCP) [13, [II)- After the conversion, each of 
the parties (nodes) knows obviously which of the 
bonds are now perfectly entangled. Using classi- 
cal communications, the parties establish whether 
there exists an infinite percolating cluster (or one 
spanning the lattice ) and who belongs to it. Then, 
by performing a series of entanglement swappings, 
it is possible to propagate entanglement between 
any two (widely) separated nodes that belong to 
the percolating cluster. We call this scheme "classi- 
cal" because it essentially maps the problem onto a 
classical bond percolation problem 19], and its suc- 
cess or failure is equivalent to the success or failure 
of bond percolation in the same lattice. Namely, 
if p > Pc, where Pc is the lattice-dependent criti- 
cal percolation threshold (or, in other words, criti- 
cal open-bond density), then entanglement between 
any two remote nodes can be established with prob- 
ability P > 0, which asymptotically docs not de- 
pend on the distance. 



surements to transform one percolation problem to 
another one. In these protocols CEP is also used, 
but is preceded by application of certain LOCC, 
which remove and replace bonds in the network lat- 
tice, resulting in a new lattice geometry with funda- 
mentally different long-range properties. The pre- 
processing in QEP may greatly improve the possi- 
bility of entanglement percolation, either by reduc- 
ing Pc, or by increasing P. 

In this paper we present a more systematic and thor- 
ough study of QEPs based on lattice transformations. In 
particular, we explore the entire parameter space with 
the addition of Monte Carlo and series-expansion meth- 
ods. Although general principles remain to be found, 
we do discover certain rules and regularities governing 
such strategies (see Sec. IIII|) . The paper is organized as 
follows. In Sec.|TT]we formulate and describe the mod- 
els and fix the notation. Section III Al discusses quantum 
networks and percolation, whereas Sec. IIIBI deals with 
LOCCs used for lattice transformations. Our main re- 
sults are presented in Sec. lIIIl which contains descriptions 
of five examples of transformations enhancing CEP: i) 
the transformation of the kagome lattice to the square 
lattice, ii) the transformation of the double-bond hon- 
eycomb (hexagonal) lattice to the triangular lattice, iii) 
the transformation of the square lattice to two decoupled 
copies of the square lattice, iv) the transformation of the 
bowtie lattice to the decoupled triangular and square lat- 
tices, and v) the transformation of a triangular lattice 
with two different degrees of entanglement into the de- 
coupled square and triangular lattices. The cases ii), iii), 
and v) were discussed already in Ref. |15|, but we present 
here more general and stronger results. We conclude in 
SecIVl 



Over the past few decades, classical percolation theory 
has seen the development of a number of quite sophis- 
ticated and powerful methods, which we have adapted 
to the questions at hand. In order to demonstrate the 
supremacy of QEP over CEP in each of the cases i)-v) 
we employ these methods combined with the methods of 
quantum information theory. An overview of these meth- 
ods is presented in the appendices: Appendix El presents 
the necessary facts from majorization theory and singlet 
conversion protocols [181] , Appendix [B] explains some de- 
tails of our Monte Carlo (MC) simulations, and finally 
Appendix [Cl deals with series expansions. 



• At the same time that we introduced CEP, we pre- 
sented several schemes that went beyond the sim- 
ple application of singlet conversion everywhere fol- 
lowed by entanglement swapping along a path. We 
call these quantum entanglement percolation (QEP) 
protocols because they use some kind of quan- 
tum pre-processing — for instance quantum mea- 



II. FORMULATION OF THE MODELS 

A. Quantum networks and percolation 

In this section, we describe the classes of quantum net- 
works and communication protocols that we investigated. 



1. Networks of bipartite states 

The networks we consider consist, prior to application 
of communication protocols, of a collection of qubits par- 
titioned into pairs, each pair being prepared in an iden- 
tical pure state |a) G C^ (g) C^. With appropriate choice 
of bases, any such state can be written 
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where the Schmidt coefficients aojQ^i satisfy ao > ai 
and ao + ai = 1- We identify these pairs with bonds 
on a two-dimensional lattice or edges on a graph, with 
two spatially-separated qubits, one occupying each end 
of the bond. At regular positions, a small set of these 
qubits are arranged near enough to one another to allow 
measurements on any subset. Such a set of qubits then 
constitutes a vertex (or node or site) which is incident to 
each edge that contributes a qubit to the vertex. This 
defines a lattice to which we apply methods of statistical 
physics. In particular, many results on the possibility of 
long-range entanglement are described using percolation 
theory [19l |20| , and depend only on the graph structure 
of the system. 



2. Classical entanglement percolation 

Given a lattice prepared as described above, we search 
for the protocols consisting of local operations and clas- 
sical communication (LOCC) that yield the maximum 
probability of achieving entanglement between two nodes 
separated by an arbitrarily large distance. The answer 
depends on the single parameter ai characterizing the 
state I a). Even with a small palette of possible opera- 
tions, finding the globally optimal solution is a difficult 
task. Instead, we search for promising classes of pro- 
tocols. The simplest protocol consists of attempting to 
convert the state associated with each bond to a singlet 
via the "Procrustean method" of entanglement concen- 
tration, which is the optimal strategy at the level of a 
single bond [13, |T^| . This conversion succeeds with prob- 
ability p = 2ai (See Appendix]^), while failure leaves 
the pair in a state with no entanglement. In this way the 
system is described exactly by a bond percolation pro- 
cess with open-bond density p. If there is a path of open 
(maximally entangled) bonds connecting two nodes, a se- 
quence of entanglement swapping measurements, one at 
each intermediate node, is then applied in order to entan- 
gle the first and last node. This protocol is the simplest 
example of classical entanglement percolation. 



3. 



Percolation theory for CEP. 



Here we review a few fundamental ideas in percolation 
theory necessary to analyze CEP. The nodes in the lattice 
can be partitioned into sets such that each node within a 



set is connected to each other node in the set via a path 
of open bonds. Such a set of nodes is called an open clus- 
ter, or sometimes simply a cluster. The central fact of 
percolation theory is that percolation processes on most 
commonly-studied lattices in dimension 2 and higher ex- 
hibit a continuous phase transition as the bond density 
passes through a critical value Pc- For p > pc there exists 
with probability one a unique (for the lattices we study 
here) cluster of infinite mass (number of nodes), while 
for p < Pc all clusters are finite with probability one. It 
follows that improving an LOCC protocol to obtain a 
small change in p can have a dramatic effect on the prob- 
ability of long-range entanglement. In the supercritical 
phase, long-range entanglement is possible, while in the 
subcritical phase it is not possible. Serving as the order 
parameter is the density of the infinite cluster 9{p) which 
we define via 
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the probability that a fixed node A (say, the node at the 
origin) is in the the infinite cluster C. When referring to 
on a specific lattice, we use a symbol representing the 
lattice as a subscript. The probability that two selected 
nodes are members of the same cluster decays roughly 
exponentially in their separation distance to an asymp- 
totic value of 0^{p). (The length scale of decay is the 
correlation length £,{p).) This means, for the problem 
at hand, that in CEP the probability that information 
can be propagated between two nodes is asymptotically 
9^{p). 



B. Transformations of lattice structure 

It was proven in Refs. [ij and [l^ that CEP is not the 
optimal strategy for establishing long-distance entangle- 
ment. The demonstration is based on applying certain 
LOCC prior to CEP. All of these pre-processing LOCC 
act on pairs of qubits, and they either transform the state 
on a given bond, or they replace two adjacent pure-state 
bonds by one, in general, mixed state bond. 

There are essentially three types of generalized mea- 
surement used: 

• Singlet conversion The optimal LOCC singlet pro- 
tocol [1^1 is used in three situations. In CEP with 
single-bond lattices we apply it directly to the state 
(1), which results in the singlet conversion probabil- 
ity p = 2ai. If the protocol is successful, the bond 
that is converted to the singlet can be used for en- 
tanglement propagation (swapping), otherwise it is 
useless. In CEP with a double-bond lattice we ap- 
ply it directly to the two copies of the state (1), 
which lives in C'^(^C'^, and has Schmidt coefficients 
Q^o, y/ctocti, V^ocki, ai, ergo the singlet conversion 
probability is 



P 
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(3) 



Finally, in all of the QEP strategies that we study, 
we also apply singlet conversion to any remaining 
untouched bonds after the lattice transformation. 

• Entanglement swapping This protocol [13], illus- 
trated in Fig. [1}l, consists of performing the so- 
called Bell measurement on a pair of qubits (b and 
c) in a node, i.e. a von Neumann measurement in a 
basis of 4 maximally entangled orthonormal states 
(in the computational basis). It allows conversion 
of a pair of adjacent singlets into a singlet connect- 
ing end points (a and d) with probability 1, i.e. al- 
lows for perfect entanglement propagation in a con- 
nected cluster of singlets. At the same time, when 
applied to a pair of imperfect states (1), it produces 
a mixed state, which, amazingly, has the average 
singlet conversion probability equal to p = 2ai. 
Unfortunately, this effect cannot be iterated. When 
applied to the mixed states, entanglement swapping 
reduces the singlet conversion probability. This last 
point places a significant constraint on our choice 
of lattice transformations. 

We use entanglement swapping in both CEP and 
QEP, but the two uses have very different effects. In 
the the case of CEP, entanglement swapping is used 
to locally move entanglement between neighbor- 
ing nodes, an operation that is repeated in hopes 
of transporting entanglement over long distances 
(This is, roughly speaking, a brute- force method). 
But with QEP, before attempting to transfer en- 
tanglement, we search for a way to selectively ap- 
ply entanglement swappings to alter the geometry 
of the lattice and hence its long-range connectivity 
properties as given by percolation theory. The goal 
of this paper and future work is to enumerate the 
rich possibilities and point a way towards a general 
description of QEP. 

• Worst case entanglement Finally, in Refs. [ij and 
lla the worst case protocol was used, which maxi- 
mizes minimal entanglement over all measurement 
outcomes. This protocol consists also of Bell mea- 
surement, but the basis is computational for one 
qubit, and corresponding to eigenstates of ax for 
the second. When applied to qubits, it produces a 
mixed state with the property that for all measure- 
ment outcomes, the resulting pure states have the 
same singlet conversion probability. We will not 
use this protocol here. 

In order to describe the lattice transformations in- 
volved in QEP, we need a more general formulation of 
the lattice than the one given in Sec. Ill A II for CEP. It 
is useful to define the percolation processes a bit more 
precisely. We begin with a graph, that is, a set of edges 
(bonds) E and a set of vertices (nodes) V. We consider 
embedding the graph in K^ in order to treat geomet- 
ric properties. In fact, the important properties don't 
change if we force the vertices to occupy points in Z^. 




FIG. 2: (Color online) a) Mapping a quantum network to a 
percolation problem. Qubits are represented by small circles. 
Each vertex (node) contains six qubits. b) using entanglement 
swapping to transform the lattice structure. The two red 
(grey) bonds in a) are replaced by the red bond in b). This 
process can be continued to produce a double-bond hexagonal 
lattice. 



The configuration of open and closed bonds can be de- 
scribed by a probability sample space ft = HeeBi^' -'^J'' 
with points cj = (tj(e) : e G i?), where u){e) takes the 
values and 1. We allow each bond to be open with 
a different probability, i.e. p = {p{e) : e ^ E). The 
appropriate measure is a product measure on f2 with 
marginal probabilities defined by fi^i^ie) = 1) = p(e), 
/ie(w(e) = 0) = 1 ~ p{e). To allow all transformations 
possible via LOCC, we take the graph to be complete — 
that is, all possible edges e = {v,w) with v,w & V are 
present (In fact, sometimes we need double bonds, as 
well.) For percolation properties, any graph with fewer 
edges can be identified with this complete graph by set- 
ting the appropriate p{e) to zero. 

Now we describe the correspondence of entanglement 
on the physical system of qubits to this percolation for- 
mulation. We denote the set of the indices of all qubits 
by A. Each qubit a G Ais assigned to a vertex v = V{a), 
with, in general, multiple qubits assigned to each vertex. 
For example, a portion of a network that can be described 
by a triangular lattice is shown in figure Fig. [2^, where 
vertex vi contains the six circled qubits. For every pair of 
qubits a,b G A with V{a) ^ 1^(&), we denote the reduced 
state by pab-, and by 5 (a, 6) the singlet conversion proba- 
bihty (SCP), that is, the probability of conversion of Pab 
to a singlet, maximized over all possible measurements. 
Then we assign p(e) ~ S{a, b) for every e = {V{a), V{b)). 
Initially, all reduced bipartite states have either p{e) = 
for separable states, or p(e) = p = 2a\ for the prepared, 
pure, partially entangled states. In Fig. [5^ the edges with 
p{e) = p are shown in black, while edges with p{e) = 
are simply absent in the diagram. Two bonds are col- 
ored red (grey) only to show that they will be replaced 
by another bond (in the sense of altering the SCP) via 
entanglement swapping. For QEP, a successful prepro- 
cessing, entanglement swapping operation such as that 
depicted in Fig. [l] alters the percolation process by set- 
ting p{{V{a),V{b))) = and p{{V{c),V{d))) = and 
p{{V{a),V{d))) — p. Likewise, in the example in Fig. [21 
an entanglement swapping measurement on qubits qi and 



92 sets pi(V{qi), yfe))) = 0, piiViq^), y((Z4))) = 0, and 
^((^(93)1 ^(94))) = P- Each of the examples discussed 
below conforms to this description. 



III. 



EXAMPLES 



We consider five examples of lattice transformations 
in this paper, each exhibiting a different combination of 
features, with implications for the analysis of protocols. 
In particular, two transformations convert one lattice to 
another, while the other three convert a lattice into two 
decoupled lattices. Four transformations involve lattices 
with single bonds, while the fifth involves double bonds. 
Quantities arising in the analysis, such as 0{p), are ana- 
lytic about p = I. In fact, the lowest-order term in the 
expansion about p = 1 is typically much larger than the 
remaining terms even relatively far from p — I, with the 
result that most of the interesting crossover behavior in 
comparing protocols occurs for smaller values of p. But 
this behavior does not appear when using the distillation 
procedure for double bonds that leads to ^ . The reason 
is that distillation produces a saturation point for 6{p) 
smaller than p = 1. This results in a much stronger dif- 
ference between classical and quantum protocols in the 
high-density regime than does singlet conversion on single 
bonds (two-qubit pure states). Four of the transforma- 
tions produce a smaller critical density on at least one of 
the resulting lattices, which gives the most pronounced 
advantage in the regime near the critical density. The 
fifth example shows that on some lattices where a partic- 
ular QEP strategy is not advantageous, allowing bonds 
of different strengths p and p' can produce regions in the 
phase space {p,p') where QEP is indeed advantageous. 
Finally, consider comparing CEP in which singlet con- 
version is applied to each bond separately, with QEP 
that results in a single transformed lattice. In every such 
case we find that QEP is better than CEP over the entire 
(non-trivial) range of p. It is an open question whether 
this is generic behavior. 

We analyze the results of the transformation of each 
lattice in three regimes: near the critical density (or den- 
sities); nearp = 1; and between these two regimes. Argu- 
ments near the critical densities typically rely simply on 
the fact that long-range entanglement is impossible on a 
lattice with density below the critical density. These re- 
sults are the most insensitive to details of the definitions 
of entanglement and connectivity. Near p = 1 we com- 
pute high-density expansions for 9(p) and related quan- 
tities (see SecO) Often the difference between the per- 
formance of CEP and QEP in this regime is marginal. It 
is important nonetheless to carry out the analysis if we 
hope to make general statements about transformations 
that hold for all p. Between the critical regime and high- 
density regime there are some techniques widely used in 
percolation theory that may be useful, such as Russo's 
formula, which is used to prove inequalities in the rate of 
change of 9{p). However, we leave these techniques for 
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FIG. 3: One of three transformations of the kagome lattice 
to the square lattice. Pairs of qubits that are subjected to 
entanglement swapping are marked with loops. 



future work and use Monte Carlo computations in the 
present paper. 



A. Kagome lattice to square lattice 

In our first example we compare CEP on the kagome 
lattice to quantum entanglement percolation consisting 
of transformation of the kagome lattice to the square 
lattice via entanglement swapping at nodes specified in 
Fig. [3l Although we do not treat them here, there are 
at least two more ways to transform the kagome lattice 
to the square lattice using the same kind of entangle- 
ment swapping. Rigorous bounds have been obtained for 
jOc (kagome) jl9| (as well as a high-precision Monte Carlo 
estimate [21|) while Pc{D) is known exactly, proving that 
the transformation gives an advantage for p lying between 
Pc(n) and the lower bound for pc(kagome). As shown in 
Table HIl the series for 9\j and 0kag are the same to the 
first non-trivial order in g = 1 — p. But the next term 
shows that 6'n > 6'kag for p close enough to 1. The MC 
data provides strong evidence that 6*0 > 6'kag everywhere 
except near p = 1, where the statistical error is too large 
to distinguish the curves. But as is evident from the 
lower plot in Fig. [3l the terms in the expansions that we 
computed dominate 0a — 6'kag already at values of p for 
which the MC data is still accurate. Thus we find that 
QEP is advantageous over the entire range of p. 



Double-bond hexagonal lattice to triangular 
lattice 



The hexagonal lattice with double bonds, each in the 
state specified by (IT|) (see Fig. [S^), was discussed in 
Refs. [ij and [l^, in which it was shown that a trans- 
formation of the lattice to a triangular lattice offers an 
advantage over CEP protocols for values of p between 
Pc(A) and another critical value defined below. Here we 
extend the analysis to determine over the entire range of 
p which of three protocols gives the highest probability of 
long-range entanglement. The first of the three protocols, 
which we call CEP I, consists of performing an optimum 
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FIG. 4: Upper plot; Monte Carlo of l-6l(p) v.s. p = 2ai. The 
solid line is on the square lattice. The dashed line is on the 
kagome lattice. Lower plot: Normalized difference between 
9q and Skag- The solid line is computed from MC data. The 
dashed line is computed from expansions from Table HIl 







FIG. 5: (a) Double-bond hexagonal lattice, (b) transformed 
into the triangular lattice (b). Pairs of qubits that are sub- 
jected to entanglement swapping are marked with loops. 



singlet conversion to each bond separately so that the 
probability of getting at least one singlet connecting two 
nodes is p' = 1 — (1 — p)^. Communication on the result- 
ing lattice is then determined by bond percolation on the 
hexagonal lattice with bond density p'. In particular, the 
critical density (see Table H is p = Pc(CEP I) w 0.4107. 
In the second protocol, we perform a more efficient con- 
version, namely distillation, on all four qubits in a double 



bond. This protocol succeeds in producing a maximally 
entangled pair with probability p" = min{ 1,2(1 — cKq)} 
(See dSl) and Appendix [XI) This results in bond percola- 
tion on the hexagonal lattice with density 



p" = min<^ 1,2 
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(4) 



We refer to this method as CEP II, with critical density 
p = pc(CEP II) « 0.358. In the third method, the en- 
tanglement swapping that maximizes SCP is applied at 
every other node (see Fig. [5]), converting the double-bond 
hexagonal lattice to the triangular lattice with bond den- 
sity p, with Pc(A) « 0.347. 

The three protocols are compared in Fig. \6\ We first 
observe that each method fails below its corresponding 
critical density, and these are known exactly. To compare 
the methods away from the critical points, we computed 
series expansions of d{p) about q = with results listed 
in Table [Til Defining d{q) = 9{1 — q) we see that 



OcEP l{q) = 0o{q^) - 1 - g' - 3q« 



(5) 



Comparing ([5]) to ^a (p) from Table |TT] we see that for 
p — + 1 the conversion to the triangular lattice is better 
than CEP I. As is evident from Figs. [6^ and [6)3, the two 
curves have the same leading behavior as p — > 1, with 
the result that the the Monte Carlo cannot distinguish 
them in this region. However, the series expansion lies 
well within the statistical error of the Monte Carlo points, 
even relatively far from p — 1 , so we can be confident that 
^A > ^CEP I even in the region where the MC cannot 
distinguish the curves. Finally, the MC clearly shows 
that ^A > 6'cEP I for smaller p, where the expansion fails. 
Turning now to CEP II, we know that 0a > ^cep ii for 
Pc{A) < p < Pc(CEP II) and from (gl) that ^cep ii = 1 
for p > 2 — \/2. It follows that there must be a crossover 
point. The MC data suggests that this occurs for p w 
0.375. In summary, we see that QEP gives an advantage 
over CEP I for all p, but that CEP II, which is more 
efficient in its use of double bonds, is better than QEP 
at high densities. 



C. Square lattice to two decoupled copies of the 
square lattice 

This transformation replaces the square lattice with 
two decoupled copies of the square lattice. To effect the 
transformation, at selected nodes the two horizontally 
opposing bonds are joined into one bond and likewise 
with the vertically opposing bonds. This procedure is 
applied at every other node, while staggering by one node 
when shifting rows as shown in Fig. [71 Because the 
transformation splits the original lattice into two disjoint 
lattices, each taking half the surviving nodes, we cannot 
compare connection between a single pair of nodes before 
and after the transformation. Rather, we consider on the 
original lattice connection between, on the one hand, a 
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FIG. 7: Doubling the square lattice. Large circles represent 
nodes at which entanglement swapping is performed. 
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FIG. 8: Monte Carlo and series expansions for doubling the 
square lattice. Lower points: 1 — Pdoub v-s. p. Upper points: 
1 — ttq^ v.s. p. Solid lines are series expansions. 



complication of splitting, this example is subtle because 
the critical density is the same on all three lattices. In 
Ref. [l^ it was shown that just above Pc{0) = 1/2, the 
transformation gives an advantage over CEP as measured 
by the choice of nodes A, A'. 

Here we examine the connectivity for all p > pc{0). 
The probability that at least one of the pairs on the de- 
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FIG. 6: Monte Carlo and series expansions about p = 1 of 
1 — 0{p) v.s. p — 2qi for various strategies on the double- 
bond hexagonal lattice, (a) Solid lines are Monte Carlo data. 
Dashed lines are series expansions, (b) Solid lines are MC. Se- 
ries expansions are indistinguishable from MC over the entire 
plot, (c) Monte Carlo data. 



pair of nodes A, A' separated by a small distance, and 
on the other hand, the same pair translated a distance 
much larger than the correlation length to a pair B, B' . 
We choose A and A' so that each one goes to a different 
lattice in the transformation. After the transformation, 
there is the possibility of connection between A and B 
on one lattice and between A' and B' on the other. We 
choose the nodes as shown in Fig. [T] In addition to the 
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FIG. 9: Monte Carlo for doubling the square lattice. Upper 



points: 1 — ^q v.s. p. Lower points: (2 



1 v.s. p. 



coupled lattices is connected is 



Pdoubic = 1 - (1 - ^n)' = ^0(2 - ^n) 



We compare Pdoubic to the probability on the original 
lattice that at least one of A and A' is connected to at 
least one of B and B' . This probability is ttq where 



TTQ = P[AeC or A' eC]. 



(6) 



(In this paper 'or' in 'X or y does not mean exclusive 
or.) We prove that the doubling is advantageous in the 
limit p ^ 1 using the power series about q ^ for 9^ 
and ttq given in Table |TT1 This follows from comparing 



TT^ = 1 - i 



j« - 36gi" 



and 



](2-0^) = l-V-32gi 



Moreover, Monte Carlo evidence suggests that the trans- 
formation improves on CEP for all p > Pc{0). In Fig. [8] 
one sees from the accuracy of the series expansion in the 
region of high-quality MC data that that the advantage 
is maintained even when the Monte Carlo becomes noisy 
near p = 1. The same MC data is not useful near Pc{0) 
for two reasons. Firstly, both Pdoub and tt^ vanish at 
the same critical point with infinite slope. Secondly, the 
systematic error due to finite lattice size and large fluc- 
tuations in cluster statistics further complicate distin- 
guishing the curves. We instead make use of an alter- 
nate expression for tt, that is ttq = 0u{2 — to) where 
u! = P[A E C\A' E C]. The condition for advantage over 
CEP then becomes 



(2-w)2 <2- 



q2 



(7) 



The MC data for u! was generated by considering the 
largest cluster in the finite lattice to represent the infinite 
cluster C even if it is not a spanning cluster. The MC 
data supports ([7]), and furthermore, shows no evidence 
of non-analyticity at pc{0), as is evident in Fig. [51 

The foregoing analysis of the doubled square lattice 
was based on the choice of A, A' shown in Fig. \7\ But 
this is not the only reasonable choice. In fact the question 
of whether the doubling transformation is better than 
CEP depends somewhat on the details of how the nodes 
are chosen. Although this is an extreme example, sim- 
ilar questions arise in analyzing other transformations. 
Thus, the ambiguity in the measures comparing the var- 
ious protocols must be addressed in future work. 



D. Bowtie lattice to decoupled triangular and 
square lattices 

This is our second example that transforms a lat- 
tice into two decoupled lattices. Figure [TO] shows the 
measurements that decouple the bowtie lattice into the 
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FIG. 10: Transformation of the bowtie lattice to decoupled 
square and triangular lattices. Loops marking pairs of qubits 
represent swapping measurements. 



square lattice and the triangular lattice. For p satisfying 
Pc(A) < p < pc(ixi) (the exact critical values are listed in 
Table HI, the transformation is obviously advantageous. 

In carrying out further analysis, the bowtie lattice 
presents a complication not present in other lattices stud- 
ied in this paper: While the other lattices are regular 
in the sense that each vertex has the same environment 
up to rotations and reflections, the bowtie lattice has 
two kinds of vertices in this sense. Notice first that all 
the surviving nodes on the square and triangular lattices 
shown in Fig. [TO] were generated from nodes of coordina- 
tion number z = 6 on the bowtie lattice, with all nodes 
of z = 4 disappearing. Each node of z = 6 on the bowtie 
lattice has four nearest neighbors of z — 6, discounting 
nodes of z = 4 (for example node A in Fig. [TO]) Further- 
more, two of these nearest neighbors are connected by a 
diagonal bond (node A") and two are not (node A'.) 

In treating the example of doubling the square lattice 
in Sec. IIII CI the fact that the lattice decouples forced 
us to consider connections between two widely separated 
pairs of nodes. We treat the present example in the same 
way, except that the more complicated local structure of 
the bowtie lattice forces us to consider a cluster of three 
nodes A, A', A" rather than a pair. In Fig. [TOb we see 
that A', A" are sent to the square lattice and A is sent 
to the triangular lattice. We consider a distant cluster of 
nodes B,B',B" related to A, A', A" by translation, and 
examine the probability that a connection exists between 
at least one of A, A' , A" and at least one of B, B', B" . On 



the bowtie lattice, this probability is ttj^ where 
TT^ = P^[A G C or A' e C or A" E C]. 



(8) 



O 




FIG. 11: Monte Carlo estimates of ttx and -PdoubA-D- The 
inset is the same data on a linear scale. Solid lines are high- 
density expansions. 



Lattice 



Pc for bond percolation 



triangular 


2sin(7r/18) f» 0.347 


bowtie 


^ 0.4045 


square 


0.5 


kagome 


fa 0.5244053 MC estimate 


hexagonal 


l-2sin(7r/18) ^0.653 



TABLE I: Pc for bond percolation on some lattices, pc(ixi) is 
the unique root of 1 — p — Qp^ + 6p^ — p ^ . All critical densities 
are exact [19l| except for Pr(kagomel[2l[|. 



Onip) 

7rn(p) 

7rnb(p) 

0a{p) 

0o{p) 

^kag(p) 



1 - g4 - 4<?'' 
1 - 4g® - ISq^ 
1- 
1- 



q' - 8q» 

6gi« + 6 



j" - 6q'^ - 21g" + 42gi5 



1 - ?=* - 3g^ - Gq'' - 25q^ 



q-" - eg'' 

,14 



4g^ 



TABLE II: A few terms in series expansions about g = of tt 
and 9 for various lattices. Oq,6a,0o, ^kag are defined via ((2)1. 
TTQb is defined via ([S}, vrg via (O , and ttcxj via ^. 





FIG. 12: a) The asymmetric triangular lattice. Solid bonds 
represent the state \a} (density p). Dashed bonds represent 
the state \a') ( density p'). b) Entanglement swapping is 
performed on pairs of solid bonds resulting in a triangular 
lattice with bond density p. The remaining dashed bonds 
form a kagome lattice which is then transformed into a square 
lattice with bond density p' as in Sec. IIII Al 



If instead we decouple the lattices via the lattice trans- 
formation, this probability is given by 



P. 



doubA-D — 



q2 



■TTDb 



2 /)2 ^ 2 



where 



TTnh = Pa[A' (^ C or A" e C]. 



(9) 



A comparison of Monte Carlo estimates and high-density 
expansions of these quantities is shown in Fig. [TT] The 
MC data shows that for small densities, the lattice trans- 
formation gives an advantage, while for higher densities 
CEP is the better protocol. The cross-over occurs at 
p « 0.425. As with the other examples, the curves repre- 
senting the series expansions suggests that no cross-over 
occurs in the high-density region. 



E. Asymmetric triangular lattice 

With this example we demonstrate that QEP can suc- 
ceed when the initial bonds are not all in the same state, 
but fail if they are in the same state. Consider the 
lattice in Fig. [T2b composed of two different kinds of 
bonds, each of the form given by ([T|). The solid and 
dashed bonds represent the states \a) (with p — 2ai) 
and \a') (with p' = 2a'i) respectively. In general, p ^ p' . 



We consider two entanglement distribution protocols: a) 
classical entanglement percolation on the original lattice 
and b) QEP consisting of creating two decoupled lat- 
tices (the square and triangular lattices) via entangle- 
ment swapping, followed by CEP on each of the resulting 
lattices. The transformation is described in Fig. [T2l Note 
that, for p' — 0, the initial lattice is a triangular lattice 
with serial double bonds and therefore the critical point 
p = •\/pc(A) sa 0.589. On the other hand for p = and 
p' 7^ the initial lattice is the kagome lattice. 

We first examine the most robust measure of a pro- 
tocol's effectiveness — the binary measure that tells 
whether long-range entanglement is possible or not. The 
phase diagrams for this example before and after the 
transformation are shown in Fig. [13l We see that there 
are regions in the phase space for which QEP is better 
than CEP, and vice versa. 

In the region in which both protocols allow long-range 
entanglement, a more detailed measure similar to those 
in previous sections is necessary. For example, for p = p\ 
we performed an analysis similar to the one used for dou- 
bling the square lattice. The two critical boundaries in- 
tersect at p = p' = Pc(A), but comparing connectivity 
just above the critical point shows that the transforma- 
tion does not improve the probability of long-range en- 
tanglement. Likewise for p = p' and p near 1, series 
expansions showed that the transformation is not an im- 
provement on CEP. We have not yet determined whether 
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FIG. 13: Phase diagram for the asymmetric triangular lattice. 
The heavy solid curve (obtained by Monte Carlo) separates 
the supercritical and subcritical regions on the original asym- 
metric triangular lattice (CEP). The dashed curve composed 
of two line segments separates the supercritical and subcritical 
regions after the transformation (QEP). After transformation, 
the supercritical region is defined to be the region for which 
at least one of the disjoint lattices is supercritical. In the light 
gray regions, long-range entanglement is possible with QEP, 
but not with CEP. In the dark gray region it is possible with 
CEP, but not with QEP. 
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APPENDIX A: MAJORIZATION THEORY 

Majorization theory has been applied to questions 
of transforming one bipartite pure state to another by 
LOCC. In particular, a theorem due to Vidal[i^| gives the 
probability that such a transformation can be achieved 
via an optimal protocol (without specifying that proto- 
col.) We state the result and apply it to the distillation 
procedure used in this paper. We must first introduce 
a certain partial order on vectors. Consider two real, d- 
dimensional vectors r and s. We define the vector r^ 
by reordering the elements of r into non-decreasing or- 
der, and likewise with s. We say that r is submajorized 
by s, denoted by r -(™ s, if 



QEP in this scenario succeeds for some other p = p', but 
this seems unlikely. 



IV. CONCLUSIONS 

In this paper we have considered the problem of es- 
tablishing long-distance entanglement in quantum pure- 
state networks on regular 2D lattices. We have discussed 
in detail several examples of quantum entanglement per- 
colation strategies that are better than the correspond- 
ing classical strategies, i.e. those consisting of a direct 
attempt to convert bonds into singlets. Our results illus- 
trate nicely the interplay between quantum information 
theory and classical percolation theory. Despite the fact 
that we do find certain rules and regularities governing 
QEP strategies, many questions remain open. For in- 
stance, we cannot say anything about optimality of our 
QEP protocols — and most probably they are not opti- 
mal. All of our protocols involve LOCC acting on pairs 
of qubits only; can multipartite LOCC, which leads in- 
evitably to creation of multipartite entanglement, help? 
We know that for sufficiently large initial entanglement, 
perfect entanglement between remote nodes may be es- 
tablished with distance-independent probability P > 0. 
Is it also possible for arbitrarily small initial entangle- 
ment? 



fc k 



(Al) 



i=0 



i=0 



for all fc = 0, . . . , d — 1. Denoting the vector of Schmidt 
coefficients of ip by A(-0), the theorem states that jV') 
can be transformed into |0) with probability p, where p 
is the largest number on [0, 1] such that X{ip) ^"^ pX{(j)). 

A simple, relevant application is the computation of 
the probability that optimal conversion of a state |'0) G 
C* eg) C* to a singlet will succeed. Application of (|A1[) 
gives p = min{l, 2(1 — ao)}. 

Now we consider distillation, defined here as the opti- 
mal protocol for converting n pure states \ai) € C^ (8) C^ 
to n — 1 product states and one maximally entangled 
state in the Schmidt bases. This operation is useful, 
for instance, in attempting to get a single, maximally- 
entangled bond from n bonds connecting two nodes. Ex- 
plicitly, 



a. 



/^^\00)^ + ■ 



with aifi > ai^i and elements of the ordered set 
(ao.o, ai.o, • ■ ■ , ttn-i.o) non-increasing. As usual, \jk)i 
is shorthand for \j)i^o '^ \k)i,i- For j E {0,1}" de- 
fine the bijective numeration k{j) — X)iLo ^^^Ji ^^'^ 
7fe = Jkia) = ao,jo"i,ji • • ■ an-i,j,._i- Then we can write 



2"-l 

E 

j:k(j)=0 



Ik [^ \j 

m=0 



mj7n/m- 



(A2) 
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In order to apply the theorem, we need to show that this 
state can be written as a bipartite state with Schmidt 
coefhcients that can be chosen to satisfy our distillation 
condition. To this end collect all the first qubits of the 
bipartite states and all the second qubits, defining 



n-l 



n-l 



\k{j))a = (g) \J^)^,o and \kij))h = (g) |j.).,l, 
i=0 1=0 

SO that (|A2[) becomes 

n-l 2"-l 



i=0 



!«.)- E VTk\HJ))a\k{j))b. (A3) 



Note that ^Lo^Tfe = ]Xi=oi'^t,o + ai,i) = 1" = 1 so 
that jASl) defines a state in C2"(8)C2" already in Schmidt 
form. The submajorization condition is 

(72._i,72._2,...,7o)^" (0,0,. 



p p 
'2'2 



for which the only nontrivial inequality is the penulti- 
mate one 1 — 7o > p/2. Thus the maximum distillation 
probability is 

p = min {1, 2(1 - ao,oai,o • • • a„-i,o)} ■ 
Note that this result agrees with the special case in (j3]). 



APPENDIX B: 



MONTE CARLO ESTIMATES OF 

e{p) AND Tv(p) 



We computed Monte Carlo estimates of 9(p) and tt{p) 
using the Hoshen-Kopelmann [22] algorithm with modi- 
fications for efficiency 23] and the Mersenne twister ^] 
random number generator. Together with series expan- 
sions and exactly known critical densities, the quality of 
the data we obtained is more than sufficient to determine 
which transformations are advantageous. We estimated 
0(p) by computing the mean density of the largest cluster 
on an L X L lattice. Near the critical density we typi- 
cally used lattices of size L = 2-5 x 10^. We took the 
value of L at the inflection point of plots of 9{p) v.s. L 
at fixed p as estimates of the correlation length ^(p). We 
typically found that £_{p) > L for p - pc < 0.0005. We 
only expect significant systematic error in this region, but 
this poses no problem because the curves are never close 
to one another in these regions. (We reformulated the 
problem when this is the case.) We computed statistical 
error, but the error bars are at most barely visible on 
the plots, so we omitted them. Exceptions are very near 
p = Pc where fluctuations in the size of the largest cluster 
become large, and near p — 1, where finite clusters are 
rare, so collecting sufficient samples to distinguish curves 
is too expensive. We discuss the effects of these errors in 
the main body of the paper. 



We wrote a single computer code to study all the lat- 
tices. The code supports lattices with vertices that oc- 
cupy points on Z^ with bonds connecting each pair of 
nearest neighbors as well as a diagonal bond from {i, j) to 
(« -|- 1, J -|- 1) . Because the connectivity properties that we 
calculated depend only on the graph structure of the lat- 
tice, we embedded the graph of each lattice in the square 
lattice plus diagonals described above. Depending on the 
lattice being modeled, some of the bonds in the underly- 
ing lattice are closed with probability one, and some are 
open with probability one, with the vertices identified in 
the embedded lattice. The graph structures of all lattices 
appearing in this paper were modeled in this way. 



APPENDIX C: SERIES EXPANSIONS 

We follow the ideas of the perimeter method [25l. l26l. l27l| 
to compute high-density series expansions of 9{p) and 
tt{p) listed in Table HTl Here we discuss the method for 
computing the series for 9{p), but our method for com- 
puting n{p) is similar. Here, a cluster is any connected 
subgraph that contains at least one vertex. We denote 
by Sq the collection of all finite clusters that include the 
vertex at the origin, and by S the partition of Sq induced 
by equivalence under translation, in other words the col- 
lection of free clusters. Choosing an enumeration aj of 
clusters in S, the probability that a randomly selected 
site is in the infinite open cluster is easily seen to be 



9{p)^i-J2^ji^-qP<i'' 



where Sj is the number of sites and bj the number of 
bonds in cluster aj, and tj is the number of perimeter 
bonds (bonds adjacent to cluster aj). Clearly, we can find 
the series expansion in q by enumerating the clusters in 
an order that is non-decreasing in tj. For instance, on 
the hexagonal lattice, only the cluster consisting of an 
isolated site has a perimeter t of size less than or equal 
to 3 (see Fig. [H|) . so that to lowest non-trivial order 
9{p) — 1 — g^. Percolation theory is a difficult subject 
precisely because the full enumeration is difficult. The 
series for 9^ (for bond percolation) was calculated by 
machine to high order in Ref. [23. Although tables of 
cluster numbers have been published, all the others that 
we are aware of are either for site percolation or for bond 
percolation with cluster size measured by the number of 
bonds rather than sites (neither of which can be mapped 
to our problem.) 

For our results we counted a few small clusters by hand, 
which is not so difficult. In practice however, we find that 
it is also essential to categorize the clusters by symme- 
try. As an example, the clusters contributing to 9 on the 
hexagonal lattice to sixth order in q are shown in Fig. 1141 
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FIG. 14: Clusters on the hexagonal lattice contributing to 9 to 
sixth order in q. Symmetries and multiplicities m are: none- 
12, one axis-6, rotation-4, two axes-3, rotation and axis-2, 
all-1. The number of contributing clusters per figure is given 
by n = ms/2, the factor of 1/2 accounting for the fact that 
each site only supports half the rotations. The contribution 
to 9 is then -ms(l - q^q*- 12. 
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